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Abstract 

. Symmetric Monge-Kantorovich transport problems involving a cost function given by a family of 

vector fields were used by Ghoussoub-Moameni to establish polar decompositions of such vector fields 
into m-cyclically monotone maps composed with measure preserving m-involutions (m > 2). In this 
■ note, we relate these symmetric transport problems to the Brenier solutions of the Monge and Monge- 

Kantorovich problem, as well as to the Gangbo-Swi§ch solutions of their multi-marginal counterparts, 
both of which involving quadratic cost functions. 



^ '. 1 Introduction 
O 



Given Borel probability measures fii, i = 0,1,,..., m — 1 on domains Q,i C R'', and a cost fmiction c : 
f2o X f^i X X S^m-i R, the multi-marginal version of Monge's optimal transportation problem is to 
minimize: 

C(Ti, ...,r„_i) / c{xf),Ti{xQ),T2{xQ),...,Tra-i{xo))d^iQ (M) 

among all (m — l)-tuples of measurable maps (Ti,T2, ...,Tm~i), where Ti : Qq ^ il^ pushes /iq forward to 
Hi for all i = 1, m — 1. The Kantorovich formulation of the problem is to minimize: 

C'{6) := / c{xo,Xi,X2, ...,Xm-i)d9 (K) 



X: 

!=! . 

among all probability measures 9 on ^Iq x Qi x ... x Clm-i such that the canonical projection 

TTi : flo X fli X ... X ilm-l 

pushes 6 forward to Hi for all i. 

Note that for any (m — l)-tuple {Ti,T2, ...,Tm-i) such that Ti^iiQ — fn for all z = 1,2, ...,m — 1, we can 
define the measure 9 = {I,Ti,T2, ...,Tm~i)^Ho on fig x fii x ... x U,„i-i, where / : f^o — > is the identity 
map. Then 9 projects to fii for all i and C{Ti,T2, ...,T„i_i) = C{9). In other words, (K) can be interpreted 
as a relaxed version of (M). 

Standard results for fairly general cost functions c show that there exists a probability measure 9 on 
fig X ^^1 X ... X Qm-i with marginals fiii * = 0,1,,...,™ — 1, where the supremum in (K) is attained. 
The natural question here is the following: 
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Problem (1): For which cost functions c. problem (K) admits a sohition 6 that is supported on a "graph", 
that is a measure of the form 9 = {I,Ti,T2, ...,Tm-i)#l-io for a suitable family of point transformations 

{Tl,T2, ...,Tm-l)- 

Whereas the case when m = 2 is already well understood, the problem when m > 3 remains elusive since 
existence and uniqueness in (M) as well as uniqueness in (K) are still largely open for general cost functions. 
There is however one important case where this problem has been resolved by Gangbo and Swi§ch [7]. This 
is when the cost function is given by 

m— 1 m— 1 

c{xo,Xl,...,Xm-l) = ^ ^ \Xi-Xjf. (1) 
i=0 j=i+l 

Other extensions were also given by Pass [16, 17, 18, 19] and by Carlicr-Nazarct [4]. 

In this note, we are interested in the symmetric versions of the Monge-Kantorovich problem. They are of 
the following type. 

Ksym = SUp<^/ c{xo,Xi,...,Xm-l)de] 6 G'Psymi^'^,!^')} (2) 



where c is an appropriate cost function and Psym{^^" , 1-^) denotes the set of all probability measures on fi™, 
which are invariant under the cyclical permutation 

Cr{xo,Xi,...,Xm-l) = ixi,X2,...,Xm-l,Xo) 

and whose marginals are all equal to a given measure fj,. Note that one can then assume that the cost 
function c is cyclically symmetric, since one can replace it by its symmetrization 



m— 1 

and in this case, one can minimize over the set V{fl"^,n) of all probability measures on fl™ whose all 

marginals are equal to /i. 

Standard results for fairly general cost functions c show that there exists £ Vgymi^"^, m), where the supre- 
mum above is attained. The natural question here is the following: 

Problem (2): For which cost functions c, problem (Kgym) admits as a solution a probability measure 9 
of the form 6 = (7, 5, 5^, S'™~^)#/U, where S is a /^-measure preserving transformation on SI such that 
S"" = J a.e. 

Problem (2) was resolved by Ghoussoub and Moameni for m = 2 in [13] and for m > 3 in [14] in the case 
where the cost function is of the form 

c{XQ,Xl,...,Xm-l) = {ui{xo),Xl) + .... + {Um-l{xo),Xm-l), (4) 

where ?i„i_i are bounded vector fields from ft R''. Their work was in the context of establishing 

polar decompositions of vector fields in terms of monotone operators, which we will briefly describe below. 
The raison-d'etre of this paper is however to make a link between the results of Gangbo and Swi§ch dealing 
with the quadratic cost (1) but for marginals of the form = cr^M for i = 0, to — 1, and the symmetric 
Monge-Kantorovich problems considered by Ghoussoub and Moameni for the cost (4) . 

Polar decompositions 

Recall that a vector field u : Cl ^ W'' on a, domain fl in M** is said to be monotone on O if for all {x, y) in f2, 

{x-y,u{x)-u{y))>Q. (5) 
A result of E. Krauss [15] states that a map u : O — >■ R*^ is monotone if and only if 

u{x) = V2-ff(x, x) for all x G fi, (6) 
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where H is a concave-convex anti-symmetric Hamiltonian on x M^. More recently, Galichon-Ghoiissoub 
[6] extended Krauss' result to the case of m-cyclically monotone vector fields, where m is a fixed integer 
larger than 2. Recall that these are the maps u from O to that satisfy for any m + 1 points {xi)'^Q in fl 
with xo = Xm, the inequality 

m—l 

E {u{xk+i), Xk+1 - Xk) > 0. (7) 

fc=0 

For that, Galichon and Ghoussoub consider the class 'Hm{^) of all m-cyclically antisymmetric Hamiltonians 
on O™, that is the set 

■Hmi^) = {H G (7(0™; M); YlZi ^ ('^'"^ W) = for all x G O™}, 

where a is the cyclical permutation a{xo,Xi, ...,Xm-i) = {xi,X2, ■■■,Xm-i,Xo). They then show that if a 
vector field u is m-cyclically monotone, then there exists a Hamiltonian H e ■Hm(O) such that 

u{x) = VmH{x, X, x) for all x e O. (8) 

Moreover, H can be assumed to be concave in the first variable, convex in the last (to — 1) variables, though 
only m-cyclically sub-antisymmetric on O™, that is H{<j'^~^('x)) < for all x e CI™. 
It is worth comparing the above to a classical theorem of Rockafellar [20] , which yields that a single- valued 
map u from O to M'' is a maximal cyclically monotone operator (i.e., satisfies (7) for every to > 2), if and 
only if 

u{x) = Vip{x) on fl, where : M'' ^ K is a convex function. (9) 

More remarkable is the polar decomposition that Y. Brenier [3] establishes for a general non-degenerate vector 
field, and which follows from his celebrated mass transport theorem. Recall that a mapping u : O — > M'' is 
said to be non-degenerate if the inverse image u~^{N) of every zero-measure A'' C R** has also zero measure. 
Brenier proved that any non-degenerate vector field u G L°°(0, R'') can be decomposed as 

u{x) = V<p o S{x) a.e. in fl, (10) 

with : M"^ — > R being a convex function and 5 : f2 — >■ fi a measure preserving transformation. 

Later, and in the same spirit as Brenier's, Ghoussoub and Moameni established in [13] another decomposition 
for non-degenerate vector fields, which can be seen as the general version of Krauss' characterization of 
monotone operators. Assuming the boundary dQ has measure zero, they show that if u G L°°{fl, R**) is a non- 
degenerate vector field, then there exists a measure preserving transformation S : il. ^ fl such that S"^ ^ I 
(i.e., an involution), and a globally Lipschitz anti-symmetric concave-convex Hamiltonian ff : R"^ x R'' — >■ R 
such that 

u{x) = 'S/2H{x,Sx) a.e. xeSl. (11) 

In other words, up to a measure preserving involution, essentially every bounded vector field is monotone, 
where the latter correspond to when S is the identity map on O. 

More recently, Ghoussoub-Moameni [14] extended this result by showing that any family ui, ...,Um-i of 
non-degenerate bounded vector field can be represented as 

Ui{x) = Wi+iH{x,Sx,S'^x,...,S"'-'^x) a.e. X en for i = 1, ...,to - 1, (12) 

where H G 'Hm{n) and 5 is a measure preserving m-involution (i.e., 5™ = /). Moreover, H could be replaced 
by a Hamiltonian that is concave in the first variable and convex in the other (m — l)-variables, though only 
TO-cyclically sub-antisymmetric. 

The proofs of the representations (11) (when to = 2) and of (12) (when to > 3) rely on symmetric versions 
of the Mongc problem and of its multi-marginal Mongc-Kantorovich version as mentioned above. We shall 
give here another formulation that is closer to the original Monge and multi-marginal Monge-Kantorovich 
problems corresponding to a quadratic cost. 
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2 Mass transport with quadratic cost in the presence of symmetry 

Given two probability measures with finite second moment /io,Mi on R**, with X := support (/Uq) and Y := 
support(/Ui), the Wasserstein distance ^2(^0, Mi) between them is defined by the formula 

VF2(Mo,m)^ = inf |^|a;-T(a;)|2d/xo(a;); Te5(/xo,/xi)| (13) 

where >S(/io, /ii) is the class of all Borel measurable maps T : X ^ Y such that T#/io = Mi) i-^-) those which 
satisfy the change of variables formula, 

/ h{y)dfxi{y) = I h{T{x))di^oix), for every h€C{Y). (14) 

Whether the infimum describing the Wasserstein distance W2(/iO)Mi) is achieved by an optimal map T is 

a variation on the original mass transport problem of G. Monge, who inquired about finding the optimal 
way for rearranging /xq into /xi against the cost function c{x) = \x\. Our cost function here c{x) = \\x\^ is 
quadratic, and the existence, uniqueness and characterization of an optimal map that we give below, was 
established by Y. Brenier. 

Theorem 2.1 (Brenier) Assume fiQ is absolutely continuous with respect to Lcbcsgue measure, then there 
exists a unique optimal map T in <S(/X0)A*i)) where the infimum in (13) is achieved. Moreover, the map 
T : X ^ Y is one-to-one and onto, /xo o,.e., and is equal to V<p no a.e on X, for some convex function 
if-.W^ 

Moreover, the Brenier map T is the unique map (up to a.e. equivalence) of the form V<p with cp convex 
such that V<^#/io = Mi- 

Now we consider the above theorem in the presence of symmetry. 

Corollary 2.2 Let fj, be a probability measure on M*^ that is absolutely continuous with respect to Lebesgue 
m,easure, and let jj, be its image by a self-adjoint unitary transformation a on R*^ (i.e., a = a* and = I)- 
Then, there exists a convex function : M"^ — ^ M such that 

yip#fj, = jl and ip*{a{x)) = ip{x) for x gR'', (15) 

where if* is the Legendre transform of (p. 

Proof: The above theorem yields a convex function : M'' — > M such that Vip^n = p,. Recall that if (p* 
is the Legendre transform of (p, then V<p* = {V(p)~^. Hence the function ip{x) = (p*{(7{x)), which is also 
convex has a gradient Vip = a* o V</5* o a, which also maps fi onto ft. By the uniqueness property, we have 
Vtjj = V<^, which means that -up to a constant- (p*{a{x)) = (p{x) for all x G R'^. 

In the rest of this section, we try to connect the above corollary to the following polar decomposition 
established in [13]. 

Theorem 2.3 (Ghoussoub-Mocmieni) Let O be an open bounded set in such that dCl has zero Lebesgue 
measure. 

1. Ifu G L°°{il,M.'^) is a non- degenerate vector field, then there exists a measure preserving transformation 
S : ^ such that = I (i.e., an involution), and a globally Lipschitz anti- symmetric concave- 
convex Hamiltonian i? : R** x R'' — >• M such that 

u{x) = V2H{x,Sx) a.e. x&Q.. (16) 

The involution S is obtained by solving the following variational problem 

sup I j {u{x), Sx) dx; Sis a, measure preserving involution onn| . (17) 

2. Moreover, u is a monotone map if and only if there is a representation as (16) with S being the identity. 
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Lagrangians and Hamiltonians associated to monotone maps 

An important example of an involutivc transformation is the transpose cr on R'' X R'' defined by (T{x,y) = 
{y,x), since the convex functions L on R'' x R'' satisfying L*{y,x) = L(x,y) for all {x,y) G R'' x R'' are 
connected to central notions in nonlinear analysis and PDEs [11]. 

Duality theory is at the heart of this concept and it is thcirefore enlightening to describe it in the case where 
R'' is replaced by any reflexive Banach space X. Recall from [11] the notion of a vector field BL that is 
derived from a convex lower semi-continuous Lagrangian on phase space L : X x X* — >• R U {+cxd} in the 
following way: for each x G X, the -possibly empty- subset dL{x) oi X* is defined as 

dL{x) := {peX*- {p,x) e dL{x,p)}. (18) 

Here dL is the subdifferential of the convex function L on X x X*, which should not be confused with 

dL. Of particular interest are those vector fields derived from self-dual Lagrangians, i.e., those convex lower 
semi-continuous Lagrangians L on X x X* that satisfy the following duality property: 

L*{p,x) = L{x,p) for all (a;,p) e X X X*, (19) 

where here L* is the Legendre transform in both variables, i.e., 

L*{p, x) = sup{(y,p) + {x, q) - L{y, q) : {y,q)eXx X*}. 

Such Lagrangians satisfy the following basic property: 

L{x,p) - {x,p) > for every {x,p) €X xX*. (20) 

Moreover, 

L{x,p) — {x,p) = if and only if {p,x) G dL{x,p), (21) 
which means that the associated vector field at x e X is simply 

dL{x) := {p e X*;Lix,p) - {x,p) = 0}. (22) 

These so-called selfdual vector fields are natural but far reaching extcmsions of subdifferentials of convex 
lower semi-continuous functions. Indeed, the most basic selfdual Lagrangians are of the form 

L(x,p) = (p{x) + (p*{p), 

where is a convex and lower semi-continuous function on X, and (f* is its Legendre conjugate on X*, in 
which case 

dL{x) = dip{x). 

More interesting examples of self-dual Lagrangians are of the form 

L{x,p) = ip{x)+ip*{-Tx+p), 

where ^ is as above, and F : X — >■ X* is a skew adjoint operator. The corresponding selfdual vector field is 
then 

dL{x) ^Tx + dif{x). (23) 

Actually, it turned out that any maximal monotone operator A is a self-dual vector field and vice- versa [11]. 
That is, there exists a selfdual Lagrangian L such that A = dL. This fact was proved and reproved by 
several authors. See for example, R. S. Burachik and B. F. Svaiter [2], B. F. Svaiter [22]. 
This result means that self-dual Lagrangians can be seen as the potentials of maximal monotone operators, 
in the same way as the Dirichlet integral is the potential of the Laplacian operator (and more generally as 
any convex lower semi-continuous energy is a potential for its own subdifferential). Check [11] to see how 
this characterization leads to variational formulations and resolutions of most equations involving monotone 
operators. 

Proposition 2.1 Let u : O — >■ R"^ he a possibly set-valued map. The following properties are then equivalent: 
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1. u is a maximal monotone map with domain O. 

2. There exists a convex self-dual Lagrangian on M'' x R'^ such that u{x) = dL{x) for all a; € O. In other 
words, 

{p,x) G dL{x,p) if and only if p E u{x). (24) 

3. There exists a concave-convex anti- symmetric Hamiltonian H onM."^ x M."^ such that 

u{x) C d2H{x,x) for a; e O. (25) 

Sketch of proof: Assuming u is maximal monotone with domain Cl, we consider its associated Fitzpatrick 
function, that is 

N{p, x) = sup{(p, y) + {q,x- y); {y, q) G Graph(u)}. (26) 
It is known and easy to see that 

N*{x,p) > N{p,x) > {x,p) for every {x,p) G M"* x M"^, (27) 

and that 

N{p,x) = {x,p) if and only if (x,p) G Graph(u). 
Now consider the following Lagrangian on M'' x R'', which interpolates between N and N*, 

L{p,x) := inf I^^N{pi,xi) + ^N*{x2,P2) + ^\\xi - Xsf + -P2||^| , 

where the infimum is taken over all couples (xi,pi) and {x2,P2) such that 

1 1 

i^'P) = 2(^i'Pi) + 2(^2,^2). 

It can be shown (see [11], p. 92) that L is a self-dual Lagrangian on M"^ x R"^, in such a way that 

N*{x,p) > L{p,x) > N{p,x) > {x,p) for every {x,p) G M'' x M'^, (28) 

which means that if {x,p) G Q x H'^, then L{p,x) = {x,p) if and only if (x,p) G Graph(w), that is when 
p G u{x). 

To show that 2) implies 3), it suffices to take the Legendre transform of L with respect to the second variable, 
i.e., 

KL{x,y) = sup{{y,p) - L{x,p);p G R''}. 
It is clearly concave-convex. The selfduality of L yields that is (at least) sub-antisymmetric, i.e., 

Kl{x, y) < -Kiiy, x) for (x, y) G R'^ x R''. 

Note that KL{x,y) < HL{x,y) = ^{KL{x,y) — KL{y,x)) is anti-symmetric, and that u{x) G d2HL{x,x) in 
Cl as in the theorem of Krauss mentioned above. 

Finally, assuming 3) that is u{x) G d2H{x,x), where H is convex in the second variable, we have for any 
x,y gCI, any p G 92-ffz,(a;, x) and q G 92-ffi,(y, y) 

H{x, y) > H{x, x) -\-{p,y- x) and H{y, x) > H{y, y) + {q,x-y). 

Since H is anti-symmetric (hence H{x,x) = 0), this implies -by adding the two inequalities- that 

> {u{x)-u{y),y-x), 

which means that u is monotone. □ 

Now we note that Monge transport problems provide a natural way to construct selfdual Lagrangians, hence 
general monotone operators. 



Corollary 2.4 Let /i be a probability measure on the product space f7 := r^i X ^2 C R'' X W^, and let p, be 
the probability measure on fl := 0,2 x fli obtained as the image of fi by the transformation a{x,y) = {y,x). 
Then, there exists a self-dual Lagrangian L : x ^ R such that VZ/#/U = jl. 
Moreover, the Monge-Kantorovich problem 

sup{/ [{x,y') + {y,x')]d9ix,y,y',xy,9 eV{n xn), proj^^ = /x, proja^ = A} (29) 

has a solution 9 that is supported on the self-dual Lagrangian manifold 

{{{x, y), {x', y')) G R^'' x M^^. y) + y') = y') + ^')}. (30) 
Proof: This is a direct application of the above corollary. The map VL then solves the Monge problem 

/ \{x,y)-yL{x,y)\^df,ix,y)=^M{[ \{x,y) - T{x,y)\^ di^{x,yy, T^fi = fl}. (31) 
Jnixn2 Jnixn2 

The self-duality of L follows from applying the above corollary to the transformation a{x, y) = [y, x). □ 
Involutions as a byproduct of niciss transport between graphs 

Suppose now that /i = (/ (g) u)^dx, where dx is normalized Lebesgue measure on a bounded domain fl in 
M'', and u : f2 — )• M'' is a vector field in such a way that /x is supported by the graph G := {{x, u{x)); x £ 0} 
and let fl = cr#/i. Now any map T pushing fj, onto fb can be parameterized by an application 5 : O — )• f2 via 
the formula: 

T : {x, u{x)) {u{Sx),Sx), (32) 
and the Monge problem between /x and p, can then be formulated as 

inf |i J [\u{Sx) — x\'^ + \u{x) — S{x)\'^] dx; S measure preserving transformation onf2| . (33) 

Assume now that -just like in the non-degenerate case- there exists a self-dual Lagrangian L such that 
VL#/[i = jl. This means that there exists S ifl ^ Cl such that 

VL{x,u{x)) = {u{Sx), Sx) for a.e. x Gfl. 

Since VL* = a* o VL o a and \7 L* (u{Sx), Sx) = (x, u{x)) a.c. wc have for a.e. x G $7, 

{x, u{x)) = VL*{u{Sx), Sx) =a*o VL{Sx, u{Sx)) = a*{u{S^x), S^x) = {S^x, u{S'^x)). 

It follows that 5 is a measure preserving involution on fi. In other words, problem (33) is equivalent to the 
problem 

sup I y (■u(x), S'x) dx; S* measure preserving involution on fij- , (34) 

which was used by Ghoussoub-Moameni to establish the polar decomposition (11). Actually, if one now 
considers the Legendre transform H of L with respect to the second variable, then as noted above, H is 
sub-anti-symmetric, can be assumed to be anti-symmetric, and satisfies, 

L{x, u{x)) -\- H{x, Sx) = {u{x), Sx) for all x Gfl. 

This then yields the polar decomposition u{x) — \72H{x, Sx) a.c., which can be then seen as a (self-dual) 
mass transport problem between the measure /x supported by the graph of u and its transpose. Unfortunately, 
the measure /i is too degenerate to fall under the framework where we have uniqueness in Brenier's theorem, 
hence the need to give a direct proof of the result as in [13]. 

In order to link the polar decomposition with the symmetric Monge-Kantorovich problem, wc note that 
measurable functions S : fl ^ fl, whose graphs {{x, Sx); x €^1} are the support of measures in 7-'^jjj(f2 x Q) 
can be characterized in the following way. 
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Lemma 2.5 Let 5 : f2 — > O be a measurable map, then the following are equivalent: 

1. The image of fj, by the map x — )• {x, Sx) belongs to 'P^y^{Q x Q). 

2. S is fjL-measure preserving and S^x = x fi-a.e. 

3. H{Sx,x) d^{x) — for every Borel measurable, bounded and antisymmetric function H on fl x fl. 

Proof. It is clear that 1) implies 3), while 2) implies 1). We now prove that 2) and 3) are equivalent. 
Assuming first that S is measure preserving such that S"^ = I a.e, then for every anti-symmetric H in 
L°°{9. X n), we have 

/ H{x,S{x))dii{x) = I H{S{x),S^{x))dii{x) = [ H{S{x),x) dn(x) = - j H{x, S{x)) dii{x), 
Jn Jn Jn Jn 

hence H{x, S{x)) d^{x) = 0. 

Conversely, if H{x, S{x)) diJ,{x) = for every anti-symmetric H, then it suffices to take H{x, y) = f{x) — 
f{y), where / is any continuous function on SI to conclude that S is necessarily /u- measure preserving. On 
the other hand, if one considers the anti-symmetric functional 

Hix,y) = \S{x)-y\-\S{y)-x\, 

then = J^H{x,S{x)) diJ,{x) = J 15^(0;) — x\diJ,(x), which clearly yields that S is an involution y^-almost 
everywhere. □ 

We now give the following variational formulation for monotone operators, which shows that they are in 
some way orthogonal to involutions. Denote by 82(^1,^1) the scit of all /x-measure preserving involutions on 
SI. It is an easy exercise to show that 52(0, /i) is a closed subset of a sphere of L'^(n,'R'^). In order to 
simplify the exposition, we shall assume that dfi is Lebesgue measure dx normalized to be a probability on 
the bounded open set fl, and /U can and will then be dropped from all notation. We shall also assume that 
the boundary of fl has measure zero. 

Proposition 2.2 Let u : O — > R'' be a vector field in Z/^(f2, R**). The following properties are then equiva- 
lent: 

1. u is monotone a.e. on fl, that is there exists a measure zero set N such that the restriction of u to 
fl\N is monotone. 

2. sup { J^{u{x), Sx — x)dx\ S £ S2{fl, m)} =0. 

3. The projection of u on S2{fl,n) is the identity map, that is 

[ \u{x)-x\'^dx = m{{ I \u{x) - S{x)fdx;S eS2{fl,ii)}. 
Jn Jn 

4. sup{J^^^{u{x),y)dTT{x,y); n G Vsymifi x fl,dx)} = J^{u{x),x) dx. 

Proof: Assume u is a monotone map and use Proposition 2.1 to write it -modulo an obvious abuse in 
notation- as u{x) = V2H{x, x) a.e. on fl, where H is anti-symmetric and convex in the second variable. We 
can write for any measurable map S, 

H{x, S{x)) > H{x, x) + {V2H{x, x), S{x) - x)). 

It follows that 

/ {u{x),x- S{x))dx= [ {V2H{x,x),x- S{x))dx> [ [H{x,x) - H{x,S{x))]dx. 
Jn Jn Jn 
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Note that H{x, S{x))dx = since S is measure preserving and S"^ = I. Similarly, J^^ H{x, x)dx = and 
it then follows that J^{u{x), x — Sx)dx > 0. Finally, by taking Sx = x, one can then see that the supremum 
in 2) is equal to zero. 

Suppose now that 2) holds. In order to show 1) i.e., that u is monotone, consider a pair xi,X2 in fi and R 
small enough so that B{xi,R) C f2 for i G {1,2}. Define the measure preserving involution Sr via 

{X — X2 + Xi if X £ B{X2 , R) 

X — xi + X2 if a; S B{xi , R) 

X otherwise. 

Since Jq{u{x),x — S}i{x))dx > 0, we have 

xi - X2, u{xi + Ry)dy- u{x2 + Ry)dy\ > 0. 

\ JBi JBi I 

Since u is Lebesgue integrable, almost every point x € SI is a Lebesgue point, which means that u{x) = 
limfl_,.o I B u(x + Ry)dy. This leads to {xi — X2,u{x\) — u{x2)) > for a.e. xi,X2 € O. 

By developing the square, it is clear that property 2) is equivalent to 3), which says that the identity map is the 
projection of u on the closed subset S2{^, l-i) of the sphere of L-^{Q, R''), that is dist(M, S2{^, A*)) = ||u — 1\\2- 
In other words, \u{x) — x\^dx = inf{ Jj^ l^(^) ~ S{x)\^dx; S € «S2(0, //)}. 

For 1) implies 4) assume u is monotone and observe that for any probability tt in ^(fi x O) with marginals 

dx, we have 

/ {u{x),y - x)d'K{x,y) = / {u{x) - u{y),y - x)d-K + I {u{y),y - x)d-K 
Jnxn JnxQ Jnxn 



< / {u{y),y - x)d'K{x,y). 
JnxQ 

Since tt is symmetric, we have then that 2 jQy;Q{u{x),y — x)dTr{x. y) < 0. The fact that the supremum is 
zero follows from simply taking the probability measure supported on the diagonal of x fi. 

Finally, note that 4) implies 2) by considering for any 5* S iS2(ri,/i) the symmetric measure dn on Q, x Cl 
that is the image of Lebesgue measure by the map x {x, Sx). Note that 

/ f{x,y)d-K{x,y) = / f{x,Sx)dx= / f{Sx,x)dx= / f{x,y)dTr{y,x). 
jQxn Jn Jn Jnxn 

a 

Back to the case of a general vector field, the above then shows that the variational problem used in the 
polar decomposition 



sup <^ J {u{x), Sx) dx; 5* measure preserving involution onfi ^ , (35) 
is nothing but a symmetric Monge-Kantorovich problem 



sup{ / {u{x),y)diT;'K eVsymi^ X ^,dx)}, (36) 
Jnxn 

where the cost function is given by c{x. y) = («(./•)• Jj)- In the case where u is monotone then the involution 
where the supremum is attained is simply the identity. 

3 Multidimensional Monge-Kantorovich Theorems 

In this section, we are interested in relating the Gangbo-Swi§ch solution of the multidimensional Monge- 
Kantorovich with quadratic cost to the following recent result of Ghoussoub-Moameni [14] . 
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Theorem 3.1 Given a probability measure fi on Cl and bounded vector fields ui,U2, ....,Um-i from fl to R'^ 
that are ji-non-degenerate, then 

1. The symmetric Monge-Kantorovich problem 

Ksym = SUp{ / [{ui{xo),Xi) + ... + {Um-l{xo),Xm-l)]d-K;Tr eVsyTn{^"',IJ')} (37) 



attains its maximum at a measure of the form = {I, S, S^, .., S™ where S is a ii-measure 

preserving transformation on Q. such that S'^ = I a.e. 

2. There exists a Hamiltonian H G 'Hm{^) such that for i = 1, ...,m — 1, 

u^{x) = V.,+iH{x, Sx, S'^x, 5^-^x). (38) 

Moreover, H could be replaced by a Hamiltonian that is concave in the first variable and convex in the 
other variables, though only m.-cyclically sub- antisymmetric. 

3. If the vector fields ui,U2, ...,Um-i are m-cyclically monotone, then (38) holds with S being the identity. 

We note first that the general multi-marginal version of the Monge-Kantorovich problem (K) where the 

probability measures /i^, i = 0,l,,...,m — 1 are given marginals on domains Q,i C M'', and where c is any 
bounded lower semi-continuous cost function c : f2o x x .... x ^m-i has the following dual problem. 

Proposition 3.1 There exists a solution 9 to the Kantorovich problem (K), as well as an m-tuple of func- 
tions (wo) ui, ..Um-i) such that for all i = 0, m — 1, 

Ui{xi) = inf (c{xo,Xi,...,Xm-i)-y^^Uj{xj)), (39) 

and which maximizes the following dual problem 



m—l „ 

^ / Ui{Xi 

i=0 •'^i 



)dpLi (D) 



among all m-tuples (uq, ui, .... Um-i) of functions Ui G L^{lJii) for which Ym^i^ Ui{xi) < c{xo, Xm-i) for 

all {xq, ...,Xm-l) ^ ^0 X ^ll X ... X r^m-l- 

Furthermore, the maximum value in (D) coincides with the minimum value in (K), and 

m — l _ 

Ui{xi) = c{xo, ..■,Xm,-i) for all {xo, ...,Xm-i) & support{9). (40) 

i=0 

In their seminal paper, Gangbo and Swi§ch [7] dealt with the case of a quadratic cost function, 

1 



c{xo,Xl,x2,...,x^_^) = J:T=o'ET=^ll\x^-XJ\^ on M'^ X M'^ X .... X (41) 



and established the following remarkable result. 

Theorem 3.2 (Gangbo-Swi§ch) Consider Borel probability measures jii on domains VLi C M.'^ , for i = 
0, 1, m—l vanishing on [d — l)-rectifiable sets and having finite second moments, and let c be a quadratic 
cost as in (4.1). Then, 

1. There exists a unique measure on JIq x fii x ... x ri„_i, where (i^) is achieved. It is of the form 
9= {To,Ti,T2, ...,Tm-i)#^J.o on Qq x fli x ... x Qm-i, where Tq = /, Tj : fio ^ andTi^fiQ = Hi. 

2. Each Ti is one-to-one iii-almost everywhere, is uniquely determined, and has the form 

T,{x) = V/*(V/i(.t)) tvhere f,{x) = |x|V2 + ip^{x) for x G (42) 

and (fi is a convex function which is related to the solutions (mOSd^ ^/ ^'^'^^ problem (D) by the 
formula ipi{x) = ^^^^^\x\^ — Ui{x) for x G W^. 
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3. Moreover, Vfo{x) = x + Tix + T2X + ... + T„_ia; for no-almost all x e R'^. 
The above result was clarified further by Agueh and Carlier [1] , who essentially established the following. 

Proposition 3.1 (Agueh-Carlier) Under the conditions of the Theorem 3.2 and with the same notation, 
we have for each i = 0, m — 1, that ^^fi is the Brenier map that pushes the measure in onto the m.easure 

_ rn—l 

V on R"^ which is the image of the optimal measure 6 by the "barycentric map" {xo, ...,Xm-i) — >■ — ^i- 

i=0 

Moreover, v is the unique minimizer of the functional v — >■ YI^Tq M^ICMij*^) where W2 is the Wasserstein 
distance. 

We now describe the situation in the case where the measures in are obtained from one measure fi by cyclic 
permutations. 

Corollary 3.3 Let ji he a probability measure on and let a be a unitary linear m-involution on M.^ , that 
is a* = and a'^{x) = x. Consider the corresponding Kantorovich problem associated to the measures 
Hi := o^^H, i = 0, TO — 1 with quadratic cost on W"^ . Then, 

1. The barycentric measure v on M.^ associated to the measure := cr^/i, i = 0, to — 1 is a-invariant. 

2. There exists a strictly convex function fo ■ — > K such that the functions fi{x) := fo{o"^~^{x)) 
satisfy 

'^fmlJ-i = i' /or i = 0,...,TO- 1. (43) 

Proof: As shown by Agueh-Carlier [1], v is the unique minimizer of the functional v X^I^o^ (Mi; ^)-: 
where fii := cr^/i. Since cr"* = 7, the uniqueness yields that v is then cr-invariant. Since now both functions 
fi and := /o o cr™~' are convex and since both V/j and VV'i = cr* o V/o o ct™"* push /Xj onto v, we get 
from the uniqueness property of Brenier maps that -modulo a constant- fi = fo° u™"'. □ 

m-cyclically monotone operators 

We shall now apply the above corollary to the case where a is the cyclic permutation cr^xo, x„i-i) = 
{xi, Xm-i,xo) on a product space X™. But first, we point to the connection with TO-cyclically monotone 
operators studied recently by Galichon-Ghoussoub [6]. They proved the following which is obviously an 
extension of Krauss' theorem to the case when to > 3. 

Theorem 3.4 (Galichon-Ghoussoub) Let ui, ...,Um-i : fl ^ M.'^ be m-cyclically monotone vector fields. 
Then, there exists a Hamiltonian H G Hm such that 

{ui{x),...,Um-i{x)) = Wx^^...^x^H{x,x,...,x) for all x G fl. (44) 

Moreover, H can be replaced by a Hamiltonian K onW^ y. {W^)"^~^, which is concave in the first variable, 

convex in the last (to — 1) variables, whose restriction to fJ™ is m,-cyclically sub-antisymmetric. The concave- 
convex function K is m-cyclically antisymmetric in the following sense: For every x = {xq, ...,Xm-i) in fl"^, 
we have 

1) + i^2,...,m(a;o ) = (45) 

m— 1 

where K2,...,m is the concavification of the function L{x) = H{a^x) with respect to the last to— 1 variables. 

i=l 

As to the connection to measure preserving TO-invohitions, they also showed that ?i : f7 M'^ is m-cyclically 
monotone /i a.e. if and only if it is in the polar set of Sm{^,l^), that is mf{J^{u{x),x — Sx)diJ,;S G 
Sm{^,lj)} = 0. Equivalently, the projection of w on ^^(fi,/^) is the identity map, i.e., 

inf{ / \u{x) — Sx\'^diJ,{x); S G Sm{'^, IJ-)} = / \u{x) — x\^diJ,{x). 
Jn Jn 
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It is easy to see that the above is also equivalent to the statement that 

sup{ / {u{xo),Xm-i)dn{x); TT eVsyni{^"',fi)} = {u{x),x)dfi{x), (46) 



and that the sup is attained at the image of fi by the map x — > (x, a;, a;), which is nothing but a particular 
case of the symmetric Monge-Kantorovich problem, when the cost function is given by c(xo, ...,Xm-i) = 
{u{xQ),Xm-i), where u is an m-cyclically monotone operator. 

Multidimensional Monge theorem on graphs and m-involutions 

We shall now make a connection between Theorem 3.1 above, which is a Monge-Kantorovich problem on 

Vsym{^"',^l) with 

c{xo,Xi,...,Xm-l) = {ui{xo),Xi) + .... + {Um-l{xo), Xm-l) 

as a cost function, and the mass transport result of Gangbo-Swicch, which corresponds to the standard 
multidimensional Monge-Kantorovich problem, i.e., with cost function 

m— 1 1 

c{yo,yi,-,ym-i) = ^ ^ Ivi-yjl'^, 

i=0 j=i+l 

and where the marginals are /Xj := ct^/Uq and fio is supported on an appropriate graph dictated by the vector 
fields ui, Um-i- 

For simphcity, we shall do this for m = 3, that is for two vector fields Ui,U2 in L°°(f2; R''). Consider 
the (degenerate) probability measure /i to be the image of Lebesgue measure dx on ft C M.'^ by the map 
P : C R'' ^ M^'' X R2d X IR2ci defined by 

X — >■ P{x) = {x, X, Ui{x), 0, 0, U2{x)). 

We now consider the 3-cycUc permutation on R^"^ x M.'^^ x K^*^ defined by 

o-((a;o,i,a;o,2), (a;i,i, a;i,2), (^2,1, a;2,2)) = ((a;i,i, xi,2), (3:2,1, X2,2), (a;o,i, 2:0,2)), 
in such a way that cr'^ = Id . 

The quadratic 3-dimensional Monge problem of Gangbo-Swi§ch applied to the measures /xq = i-ii = a^ii 
and fjb2 = CT^M becomes the problem of minimizing 

/ {\\P{x) - aP{Six)f + \\aP{Six) - a''P{S2x)\\'' + \\a''P{S2x) - P{x)f] dx 
Jn 

over all measurable maps (5*1, ^2), where Si : ft ^ fl is measure preserving for i = 1,2. 
The Kantorovich formulation of the problem is then to minimize 

C{n) := / {\\P{x) - aP{y)f + \\aP{y) - a'P{z)f + \\a'P{z) - P{x)f}dn{x,y,z) 
over all probability measures d7r(x, y, z) whose 3 marginals are Lebesgue measure. In other words. 



= / {\x-u,{y)\' + \x\' + \u,{x)\' + \u2iy)\' + \y\' + \u2{x)-y\ 



\\x - uiyy)\ -r PI -r \ui(x)\ -r \u2\y)\ T \y\ -r \u2\xj - x/'^ 

/(M<i) 

|2 I L, /„m2 , ii2 I u, „|2 I I , „, /^m2 , \„,\2 



+\My)r + \M^)r + + \My) - A' + \y- + w 

+\x\^ + \x- U2{z)\^ + \ui{x) - z\^ + \z\^ + \ux{z)\^ + \u2{x)\^}dTT{x,y,z). 

Since the integrals of x^,2/^,2:^,Mi(a;)^,ui(j/)^,ui(2;)^,U2(a;)^,U2(y)^,U2(-2)^ against the given marginals of tt 
are given constants, the above problem amounts to minimize 

/ {\x - ui{y)\^ + \u2{x) - y\^ + \u2{y) - z\^ + \y- ui{z)\'^ + \x- U2{z)\^ + \ui{x) - z\^} d'K{x,y,z), 
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or, for the same reasons, to maximize 

D{Tr) := / {(wi(y),x) + {u2{x),y) + {u2{y),z) + {ui{z),y) + {u2{z),x) + {ui{x),z)] d-K{x,y,z), 

J(R<i)3 

which is exactly the problem (Kgym) where the cost 

c{x,y,z) = {U2{x),y) + {ui{x),z) 

has been symmetrized. Consider now the optimal maps Ti , T2 obtained by Gangbo-Swicch, that is Ti = 
V/* o V/o pushes /xo to jJi = cr^/Uo and where fi is strictly convex for i = 1,2. By Corollary 3.3, we have 
that /o = /i o CT* for i = 1,2. 

Note now that there exist measure preserving transformations , ^2 on f2 such that the optimal maps 

fix := Ti{Px) = V/i* o Wfo{Px) = aP{Six) = {ui{Six),0,0,U2iSix), Six, Six) 
maps dx onto /xi, while 

Tax := T2{Px) = V/2* o V/o(Pa;) = a'^ o P{S2x) = (0, U2{S2x), S2X, S2X, ui{S2x), 0) 
maps dx onto fi2- 

Let now T2,i := V/| o V/i in such a way that T2 = T2,i o Ti, and let ^2,1 be a measure preserving on Cl such 
that 

T2,i{ui{y),0, 0, U2{y), y, y) = (0, U2{S2,iy), S2,iy, S2,iy, ui{S2,iy), 0), 

so that 

T2,i o Tix = (0, ^2(5*2,1 o Six), 5*2,1 ° Six, 5*2,1 o 5ix, wi(52,i o 5i.t), 0). 

Since T2 = T2,i o Ti, we have that 52, 1 o 5i = 52, and since V/o = o"^"* o V/i o cj' for i = 1,2 and 
^2,1 o (j{Px) = V/2 o V/i o a{Px) for a.e. x G fi, one can easily verify that 

52,1 = 5i := 5, 52 = 52,1 o 5i = 5^ and 5^ = /. 

In other words, the points (to, T'l^O: 72io) s-re such that 

to = P{x) — {x, X, ui{x), 0, 0, U2{x)), 

Tito = aP{Sx) = {ui{Sx),0,0,U2{Sx), Sx, Sx) 

and 

T2to = a'^PiS'^x) = {0,U2iS'^x),S'^x,S'^x,ui{S'^x),0) 

where 5 is a measure preserving transformation such that 5^ = /. 
The convex function $o(i) = /o(t) — sl^^l is such that 

V$o(to) = Tito+T2to = {ui{Sx),U2{S'^x), S'^x, S^x + U2{Sx), Sx + uiiS'^x), Sx). 

Define now the convex Lagrangian 

L{x, y, z) = ^o{x, X, y, 0, 0, z) =: 

We then get 

VL{x,y,z) = {Do,i<^it) + Do,2^t),Di,i^{t),D2,2Ht)) 

and 

VL{x,ui{x),U2{x)) = {ui{Sx) + U2{S'^x), S'^x,Sx). 
Let now H be the Legendre transform of L with respect to the last two variables, that is 

H{x, y, z) = sup{{y,p) + {z, q) - L{x,p, q);p e g € M'*}. 

L is clearly concave in the first variable, convex in the last two variables and 

L{x, ui{x), U2{x)) + H{x, S'^x, Sx) — {ui{x), 5^x) + (u2{x), Sx), 
and in other words, {ui{x),U2{x)) = 'V2,3H{x, S^x, Sx) for all x £ fi. 
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